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We re-examine the re-scattering mechanism for the X(3872), as a candidate for the 2P charmonium 
state Xci(2P), decaying to J '/t/jp(w) through exchanging D^*' mesons between intermediate states 
D(D) and D* (D* ) . We evaluate the dispersive part, as well as the absorptive one, of the re-scattering 
amplitude and find that the contribution from the dispersive part is dominant even when X(3872) 
lies above the threshold of the neutral channel th n = m D o + m D ->o. We predict R p / U — 1 for the 
mx region scanned by experiments. Meanwhile, we also estimate the rate of X — > D°D ir°. Our 
results favor a charmonium interpretation of X(3872) when it lies slightly below the threshold of 
D°D*°. Furthermore, we evaluated the width of X — » J/ipp with the help of a phenomenological 
effective coupling constant gx, and find the total width of X(3872) to be in the range of 1-2 MeV. 
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I. INTRODUCTION 



B(X -^jJ/ip)/B(X 



-J/ip) = 0.14 ±0.05, (4) 



In recent years there have been a number of exciting 
discoveries of new hadron states (for a recent review, see, 
e.g. [![). These discoveries are enriching and also chal- 
lenging our knowledge for the hadron spectroscopy, and 
the underlining theory for strong interactions. Among 
these new states, the X(3872) may be the most mysteri- 
ous one, which was first discovered by the Belle collab- 
oration j2[] in the invariant mass spectrum of J/-07r + 7r~ 
in the decay B + — > Jjipn + ir ~K + , and confirmed soon 
by BaBar U, CDF @] and DO @ collaborations. The 
world average mass is mx — (3871.2 ± 0.5) MeV and 
the width is T x < 2.3 MeV at 90% C.L. [6], which is 
consistent with the detector resolution. The dipion mass 
distribution in J /ipiT + TT~ seems to favor a p resonance 
for the dipion structure. This implies the C-parity of 
X(3872) is even, which is finally confirmed by the mea- 
surement of V(3872) jJNj @, @|. The angular dis- 
tribution analysis by Belle Q favors J PC = 1 ++ . Anal- 
ogous analysis fio] ] and the analysis for the dipion mass 



spectrum 



lysis 



by CDF collaboration allow J 



PC 



r 



and J PC = 2 + as well. The recent observation of the 
near threshold decay X -> D D°n (with a slightly 
higher mass of about 3875 MeV) by Belle may favor 
J PC = 1 ++ but can not rule out J PC = 2~+. More- 
over, Belle also see the sub-threshold decay X — » ujj/tp 
in B+ -> J/iPn + Tr-7r K + Q. So far, for the X(3872) 
four decay modes have been observed with following 
fractions [iH El] 



B{B- 



^X) X B(X -> TT+TT" J/lp) 



(1.14 ±0.20) x 10" 



B(X -> J/lfjTT+n-TT ) 

B{X ^ J/iPir+Tr-) 



1.0±0.4±0.3, 



B(B -> KX) x B(X -> D°D°tt°) 
= (1.22±0.31±g'jg)x 10- 4 , 



(1) 



(2) 



(3) 



where the experimental value for X — > ^yJ/ip is taken 
from the Belle measurement |7[ , while the observed value 
of about 0.25 by BaBar is somewhat larger Q. 

For convenience, we define the following ratios and 
their values can be deduced from (J5J) and ©: 



- fp 



= 1.0 ±0.5, 



R 



p/DDir 



Dtt 



= 0.10 ±0.05, 



(5) 



(6) 



where Ti denotes the width of decay X — > i with i = 
ipp, tpu and D°D°ir°, respectively. 

Because of the closeness of mx to the threshold 
M D o B , = 3871.81 ± 0.36 MeV [Hj], many authors iden- 
tify the X(3872) with a molecule of D°D*° + c.c. in 
S-wave (HI, a loosely bound state of charmed mesons. 
This is certainly a very attractive interpretation, which 
also gives a natural explanation of the J of X(3872), 
and predicts R p / U ~ 1 (see Ref. [HJ]) as well. Thus, 
the molecule becomes the most popular interpretation 
for the X(3872). However, it seems to be difficult for 
the molecule models to account for the large production 
rates of X(3872) at B-factories and the Tevatron [l6| un- 
less B(X — > J/4>p) is large [13], which, however, seems to 
be in contradiction with Furthermore, the molecule 
model predicted the decay into J/ipp to be much su- 
perior to that into D°D°tt°, but this seems not to be 
supported by the experiment. Moreover, the molecule 
model predicted that the production rate of X(3872) in 
B+ -> XK+ is much larger than that in B° -> XK° [l^, 
but the Belle data show that the rate of B° — * XK° 
with X — > D°D°ir is approximately equal to that of 
B + — > XK + though the errors for the measurements are 
large 12j. So, it might be useful to try other possible 



interpretations for the X(3872). 

Motivated by the large production rates in B decays 
and in pp collisions at the Tevatron, we suggested that 
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the AT(3872) be a J PC = 1++ x c i(2P) charmonium- 
dominated state [T^] . This possibility has also been sug- 
gested by Suzuki with detailed discussions on its de- 
cay properties (l9j . In this charmonium picture, the 
large rate of B — * X(2>%72)K , which is comparable to 
(not much less than) B — * XciiXP)K , and the sim- 
ilarity in production between X(3872) and ip(2S) ob- 
served by the CDF and DO collaborations at the Tevatron 
can be well understood Leaving the mass prob- 

lem [2(], HH [H, [23[ alone, the most difficult problem 
of this assignment is how to explain the observed large 
isospin violating effect expressed by Rp/^ ~ I, since the 
state Xci(2P) is an isospin scalar. Suzuki estimates 
that R p / U ~ 1/2 in an semi-quantitative way in which 
both J/i[>P and J/tpcu are produced through the D(D*) 
exchange between DD* pair, and the large isospin vio- 
lation can be accounted for by the mass difference be- 
tween neutral and charged DD* thresholds and the large 
difference between the phase spaces of X — ► J/tfjp and 
X — > J/tpui. Recently, the ratio R^mp-w was studied in 
a similar but more quantitative way [2J|, and was pre- 
dicted to be R p /DD7r ~ fO~ 6 -IO -4 , which is far smaller 
than the experimental data in 

The estimation of R p /ddtt given in Ref. [Hf is based 
only on the imaginary part of the amplitude A{X — > 
D°D*° + c.c. — > J/4>p). The quantity |Iim4| 2 can be un- 
derstood as the probability of finding the final state J / tpp 
through re-scattering of the real D°D*° +c.c. pair in per 
unit of final-state phase space. Then InxA is proportional 
to the phase space factor of X — ► D°D*° + c.c, which 
is small or even zero, since the mass of X is taken to be 
very close but above the D°D*° threshold. As a conse- 
quence, the value of R p /ddk given in Ref. [24[ is small. 
Furthermore, since the charged channel D + D*~ + c.c. is 
forbidden by phase space, this mechanism will predict 
almost equal amplitudes for J/ipp and J/ipuj, and then 
result in a large R p / U of order 10 (or even larger) due 
to the difference between J/tfjp and J/tpw phase spaces 
(see, e.g. 0). 

On the other hand, in contrast to the imaginary part, 
the real part of the re-scattering amplitude, which repre- 
sents the effects of virtual intermediate states DD* , may 
be dominant in this case, since it does not suffer from the 
phase space suppression for producing a real DD* pair. 
Moreover, the real part may also give a moderate value 
of Rp/uj, since both the neutral channel D°D*° + c.c. and 
the charged channel D + D*~ +c.c. can contribute through 
these virtual effects. 

In this paper, we re-explore the re-scattering mecha- 
nism and evaluate the real part of the amplitude as well 
as the imaginary one. With a reasonable choice for the 
phenomenological parameters, we find that the experi- 
mental data in §5§ and ^ can be explained quite well if 
X(3872) is a 2P charmonium state x c i(2P) lying below 
the threshold of D°D*°. 
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FIG. 1: The decay diagrams for X(3872) -> D*°D° + c.c. -> 



II. THE MODEL 



A. X -> J/ipp(u>) 



In the re-scattering mechanism, the decay X — > 
J jippiui) can arise from exchange of a D^*> meson be- 
tween D(D) and D*(D*). The Feynman diagrams for 
X — > D°D*° + c.c. — > J /ipp(uj) are shown in Fig. Q] and 
those involving charged intermediate states can be eas- 
ily obtained through replacements of D°D*° by D + D*~ 
and D°D*° by D~ D* + in Fig. [TJ Since mx is very close 
to the threshold M DB , and the X(3872) couples to DD* 
in an S-wave, we assume that the re-scattering contri- 
butions are dominated by the intermediate states DD*, 
and those arising from higher exited D meson states are 
neglected. 

Assume that X(3872) is the x c i(2P) state, and all the 
vertexes in Fig. [T] are determined by the effective La- 
grangians, which are constructed based on the chiral and 
heavy quark spin s ymm etries and parity conservation (for 
a review, see Ref. |25j). These Lagrangians read [26l. l29j 
(for convenience here we use the same notations and sym- 
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bols as in Ref. [H|): 

Cx = gxX»(DD*; -D^D*), 



(7a) 



CpDD = ig^DDi>„ (d^DD^ - Dd»D^) , (7b) 
C+D'D = -g^De^d^(d a D*^+Dd a D*^, (7c) 

+^D*"df,D*J - Tl) v d ll D* v B*rt } , (7d) 

C DDV = -igDDvDtd^WYp (7e) 
C D , DV = -2f D , DV e^ a p{d^))(D^ a D*^ 

-D^Td*), (7f) 
£ D , D , y = +ig D *n*vD* v ^ tl D*J{W 1 )) 

+±if D * D *vD*l{d»N v - 0"V*)*I#, (7g) 



where the indexes i,j in (|7el|7g[ ) represent the flavors of 

light quarks, i.e. D^*> =(ZjM° , flW", L>i* )_ ) T , and they 
are hidden in (|7all7d[) . V is the 3x3 matrix for the nonet 
vector meson, 



V = 



£ a 

V2 + V2 



K* 



P -^5 + ^ K*° 
K*~ K*° 6 



(8) 



We assume that chiral symmetry is preserved in ([7|. 
That is, the coupling constants are blind to the flavor and 
there are no isospin violations at the Lagrangian level. 
All the coupling constants will be determined in the next 
section. However, it is necessary to emphasize here that 



the determinations will not account for the off-shell effect 
of the exchanged D{D*) meson, of which the virtuality 
can not be ignored. As shown in Ref . 24| . such effects can 
be accounted for by introducing, e.g., the monopole form 
factors for off-shell vertexes. Let q denote the momentum 
transferred and rrii the mass of exchanged meson, the 
form factor can be written as [2J, [2al 



A 2 - mt 



A 2 - q 2 

and the cutoff A can be parameterized as 



A(m, 



aA 



qcd- 



(9) 



(10) 



We are now in a position to compute the diagrams 
in Fig. Q] If the X(3872) lies above the D°_D*° thresh- 
old, in the process X(p x ,ex) -* D°(pi) + D*°(p 2 ,e 2 ) -> 
J/ip(p3, £3) + p(uj){p4, £4), where the momenta p and po- 
larization vectors e are denoted explicitly for the mesons, 
we can calculate the absorptive part (imaginary part) of 
Fig. [Ha) and find it to be given by 



Abs(a„) 



\Pi\ 



dQA(X -> D°D*°) 



32Tr 2 mx 

xA a (D°D*° J/M")), (11) 
where p\ is the 3-momentum of the on-shell D° meson in 
the rest frame of X(3872) and the subindex "n" denotes 
the contribution coming from the neutral channel (for the 
charged channel we use " c" ) . Analogous expressions can 
be found for diagrams (lb-Id), and the absorptive parts 
of diagrams (le-lh) are the same as those of diagrams 
(la-Id), respectively. Explicitly, the absorptive parts of 
diagrams (la-Id) are given by 



J 



Abs(a„) = -^T^r — / dn{A^/2g x g^ DD f D * DV ) ^-^^\ ' (Pi ' ^lY^^p^V^x^ 

tj^iTT 777ji£" J y 777^ 

, \pi\ f , nr /x ^ 2 (m 2 ,q 2 ) „ »„ a 

Abs(6„) = — -J / dil{V2gxg^D'DgD'D'v) — 2 — £ ^ a l 3 P3 e 3 Pi 

tjZilT 777^ J 

x • Wx - — Z2—P2\ + 2r[p4 ■ ex -3 Je/ - 2r[e 4 • e x -3 J^2 

I "'2 m 2 ''*2 

Abs(c„) = — [ dniV^gxg^D'DgDDv) ^ \ pi ' £4)£^a/3P3 e 3>2 e X: 

327T z mx j qr* — 

A , /, n f , n/ „ , T 2 {m 2 ,q' 2 ) ii *v a 

Abs(d„) = -— — / dlt(2V2g x g^ D » Dr f D *DV) — 75 n~ e^apP^i Pi 

d2-K z mx J q — m 2 

*x /3 . r (gg ■ ex)(p2 - Pa) ] *p w * s , (pa ■ ex)(p2 • e|) i p 
x <^ 2(p 2 • e 3 )4 + [p 3 • e x J e 3 p - [{e 3 ■ e x ) + -5 JP3 



(12) 

where the ratio r = /d* ny /g_o* ny and the momentum transferred q = p% — p±, q' = p±— p\. The imaginary parts of 
the charged channel amplitudes are the same as (fT2"|) for X — > J/i/joj and of opposite signs for A — > J ftp p. The total 
absorptive part of the neutral (charged) channel amplitude can be obtained by a simple summation and read 

Abs„ (c) = 2[Abs(a n(c) ) + Abs(6„ (c) ) + Abs(c n(c) ) + Abs(d„ (c) )], (13) 
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where the factor "2" comes from the equality of contri- 
butions from diagrams (la-Id) and (le-lh). 

The amplitudes in (|12[) are almost equal to those given 
in Ref. except that some minor errors in Ref. [HI 
have been corrected. All these amplitudes are propor- 
tional to the phase space factor 



mx 



•x 



(mi + m,2) 2 )(m 2 x — (mi — m 2 ) 2 ) 



2m 2 x 



I m x — {mi + m 2 ) 



2m x 



(14) 



which is very small even if X(3872) is above the D°D*° 
threshold. 

In the case that X(3872) lies below the D°D*° thresh- 
old, the absorptive part (imaginary part) vanishes, and 
the dispersive part (real part) of the re-scattering ampli- 
tudes will play the role in the decay. 

The dispersive part of the re-scattering amplitude can 
be obtained from Abs„ and Abs c via the dispersion re- 
lation peh, m 




D'° 
(b) 



D° 



FIG. 2: The diagrams for X(3872) -> D°D°ir. 



Breit-Wigner distribution as 



i r {mx - m3)2 r° p {t)m p r p 



7T J m ^ (t-m 2 p ) 2 +m 2 p r 2 p dt > 

uj(7r + TT-7r°)J/iP) = 

i /-(™^-™3) 2 TlJt^r, 



(t-ml) 2 +mlTl 



dt, (16) 



where m 2 -n and m^^ are the experimental cutoffs on the 
(7r + 7r~) and the (7r + 7r~7r°) invariant masses respectively, 
and r p ( w ) denotes the total width of p(u>). 



B. X — D°D 7Y° 



Dis(mx) - 



2 Abs„(s') 



ds'- 



thl s 



•x 



thl s 



AbSc(5 2 /} ^ ),(15) 



•x 



where th n = m D a + m D *o and th c = m D ± + m D *^ are 
the thresholds of neutral and charged channels respec- 
tively, and the contributions arising from higher chan- 
nels are neglected in (fT5|) as we have mentioned before. 
Unlike the absorptive part, the dispersive contribution 
suffers from the large uncertainties arising from the com- 
plicated integrations in (JTSJ) . Since the absorptive parts 



Abs„( c )(s) falls off as s increases, it is reasonable to 
choose a cutoff for the integration to make a numeri- 
cal estimation. Following Ref. [l9j], we choose the cutoff 
around s max = Am 2 D „ , which can shut the widows of 
higher channels automatically. 

It is worth emphasizing again that for X — > J/ipu; 
the contributions from neutral and charged channels are 
nearly equal and share the same sign, while for X — > 
J/tpp they almost cancel each other. This is not surpris- 
ing since the explicit chiral symmetry is maintained in 
the effective Lagrangians in ((7|) . So if we neglect the ab- 
sorptive part, the isospin violation, which is mainly due 
to the difference between th n and th c and that between 
the thresholds of J /tpp and J/i/joj, seems to be too small 
to account for the experimental data in ([5]). However, the 
large difference between the phase spaces of X — » J/ipp 
and X — > J/tpu due to the large width of p resonance 
may result in a favorable prediction for R p / U . To achieve 
this, we smear the width I 1 ^ / Jt), which is obtained 
through the re-scattering amplitude in the narrow width 
approximation (NWA), over the variable t = m\ by the 



If X(3872) lies above the D°D*° threshold, the width 
of X — > D°D°Tr a can be given by 

r(X^D D Tr ) = 2r{X^D D* )Br(D* ->D°Tr ), (17) 

where the branching ratio Br(D*°— > D n°) is known Q 
and the width T(X -> D°D*°) can be easily obtained 
from C x in the NWA: 



r(A^L> u D* u ) 



9x\Pi\ 
2Ai:m 2 x 



\P1 



&irm 2 x ' 



(18) 



where the 3-momentum pi is the same as in ([lip . 

On the other hand, when mx is below the threshold 
th n , it can decay to D°D°t: through virtual D*°(D*°) as 
illustrated in Fig. [2] [27[. Here, we need another effective 
Lagrangian to describe the D°D* n° coupling [25l [27j : 



Cn. 



D'Dn 



.gD'D-rr 



Dfd^TT a D° - D°d fl Tr D* p ) . (19) 



Then the amplitude for X (px , £x) ^> D° D°ir° (kz) reads 

iV2gxgD*DTT 



iM = i(M a + M b ) = 



q 2 - m 2 D , Q + im D *oT(D*°) 



x L —2 (*8«ex)J, 



(20) 



where q = px — Pi is the momentum transferred and 
T(D*°) is the total width of D*°. It can be verified 
that the amplitude M. generates the same width as that 
given in (fT7|) in the limit mx — th n T(D* ). However, 
the validity of (|20[) and HHJ) are questionable in the near 
threshold region where \mx ~ th n \ ~ r(D*°) since the 
perturbation calculations are known to be invalid in this 
region. 
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TABLE I: Parameters used in the calculations. 



J/4' 
■ 4 


D gv 


P 


A 


9 


9i>DD 


5.9 


0.9 


0.56 Gev" 1 


0.6 


8 




g-ipD'D* 


u 


h-QCD 


a 




4.3 Gev" 1 


8 


132 MeV 


220 MeV 


4 



FIG. 3: The diagram for calculating the matrix element 
{D\crf ll c\D) in the VMD mechanism. 



III. NUMERICAL RESULTS AND 
DISCUSSIONS 

A. Parameter determinations 

Since the numerical results are indeed sensitive to some 
of the parameters introduced above, we need to explain 
how we determine these parameters. 

The coupling constants in Eq. (|7el|7g[ ) are universal for 
p and u>. They can be related to the standard parame- 
ters in the so-called heavy meson chiral Lagrangian [25| 
through the relations [2J, |26| 



9dd\ 



9d* d*v 



V2 



f D* DV 



V2 



The values of gv, (i, and A used here are the same as 
in Ref. [24[ and are listed in Tab. HI Similarly, the 
coupling constant in (|19D can be related to the well- 
known parameter g |25ll27| through the relation gjj* ^ — 
2\/niDmD'g/ /tt, where / w is the decay constant of it and 
the value of g can be determined by the measurement of 
the width of D* + As a byproduct, we can estimate 
the total width T(D*°) » 0.07 MeV by using the value 
of g listed in Tab. |T] together with the branching ratio 
Br(£>*° -► D°7r°)=(61.9 ± 2.9)% [|]. 

The coupling constant g^DD in (|7b[) can be estimated 
by the vector meson dominance (VMD) mechanism [H, 
|29| ]. Considering the matrix element {D\cq ^c\D) , it can 
be represented in Fig. [3] in the VMD mechanism, where 
the circle vertex is determined by Eq. (|Tb[) and the box 
vertex is related to the J/tp decay constant through 
the matrix element (0|c7^c| J/ip(p, e)) = f^m^e^. Then 
at the normalization point, where the initial and final 
D mesons have the same 4-velocities, we can determine 
g-ipDD = Tnip/fy ~ 8, which is consistent with the pre- 
diction of QCD sum rules [3(J. Other coupling constants 
in (fTc)) and (|7d[) can be estimated through heavy quark 
symmetry relations: g^D'D* = m D g^, D * D = g^ DD 
One should notice that use of VMD here does not mean 
that all higher resonances give contributions far smaller 
than those from J/tp, but it lies on the argument that 
these contributions tend to cancel [1^, [3l|. For exam- 
ple, the analogous effective coupling constant govern- 
ing ^(3770) — > DD decay is about 3 times larger than 
g^DD dll- 

The coupling gx is not involved in the ratios R p / U 



and Rp/DDTT- However, gx is important for deter- 
mining the decay widths and clarifying the proper- 
ties of X(3872). Assuming that X(3872) is a pure 
charmonium 2P-state x c i(2P), we parameterize gx — 
2y/2mj}omu'omxgi(2P), where gi(2P) is the coupling 
constant governing the interactions of 2P charmonium 
states with £)(*)£)(*) $2%. In Ref. [H, the IP partner of 
gi(2P) is estimated in a similar way to that for g^DD- 
One only needs to replace the vector current = cy^c 
by the scalar one S = cc, and the J ftp by the XcO i n 
Fig. [31 and the result is [32| 



giO-P) = 



XcO 



1 



3 /; 



(21) 



XcO 



where the 
(0|6c|xco(p)) 



decay 



constant f Xc0 if 
m v „ n . Using f Xa0 



defined by 

fxco m xco- Using fxao = (510 ± 40) 
MeV estimated by the sum rule analysis [32|], one can 
get ffi (IP) S3 2.1 GeV- 1 / 2 . As we have mentioned, for 
the charm systems g\ (2P) should be of the same order as 
gi(lP). Then, the value of gx can be estimated through 
gx ~ 2\/2m£)0 77j D ,omx<?i(lP) ~ 23 GeV. On the other 
hand, the effective interactions between charmonium 
and DW^W can also be estimated by the quark pair 
creation models (33[. From available calculations in 
Refs. [IS H3 together with Eq. (fTf]). we can deduce the 
effective coupling at hadronic level gx — 8-15 GeV when 
Smx = mx — th n varying from 80 MeV to 0.5 MeV. 
Based on the two estimates mentioned above, we will 
choose gx = 20 GeV in our calculations. This should be 
a reasonable choice for the coupling which describes the 
Xd(2P) decay to DD*. 

Since the virtuality of exchanged meson in Fig. [1] is 
always larger than 1 GeV 2 , the amplitudes in (fT2"]) are 
sensitive to a when a < 3. The authors of Ref. [24| 
choose a= 0.5-3.0. In Ref. (32| |. it is argued that the 
value of A in © can be around 3 GeV, which corresponds 
a ps 5. In our calculations we choose a = 4. 

For the charm meson masses we take m D a = 
1864.847 ± 0.178 MeV [3 and m D ,o ~ m D o = 142.12 ± 
0.07 MeV 0, so the threshold i/i„=3871.8MeV. For 
other mass and width parameters, we refer them to 
PDG2006 [6| . The cutoffs on dipion and tripion invariant 
masses in ((TBI are taken to be the same as in the Belle 
experiments^, 0] : 



m 27r = 400MeV, m 37I = 750MeV. 



(22) 



_^ Wx, (MeV) 

3871.0 3871.4 th„ 

(b) 




(d) 

FIG. 4: The mx-dependence of R p / U and R p /ddtt- ( a ,b) for 
the region m x = 3870.8-3871.8 MeV and (c,d) for m x = 
3871.9-3874.5 MeV. 



B. Numerical analysis 

Our numerical results for the mx-dependence of R p / U 
and Rp/DD-n are illustrated in Fig. |U In the below- 
threshold region where m x = 3870.8-3871.8 MeV, 
Eq. ([20)1 is used to deduce the width Too^. For the 
region above the threshold th n with mx = 3871.9-3874.5 
MeV, we use Eq. (JTTJ) and (HSJ) to calculate T DD7T . The 
contribution from absorptive part of the re-scattering 
amplitude is not involved in Fig. [4J since it is numeri- 
cally far smaller than that from the dispersive part even 
after the phase-space smearing in (|16[) . For compari- 
son, we also choose the naive cutoffs m2n — 2m„ and 
m^ = Sm^ to evaluate the integrations in (|16p . and the 
results are shown in Fig. 4(b) and Fig. 4(d). As usual, we 
use the central values of the parameters given in the last 
subsection. 



FIG. 5: The width of X — ► J/ipP- TJ/ip arises from the 
real part of the re-scattering amplitude and F^ p from the 
imaginary part. Here gx = 20 GeV is used. 

From Fig. BJa,c), one can see that in both regions of 

m x , 

R p/uJ = 1.0 ±0.3, (23) 

which is consistent with Eq. ([5]). This result is sensitive 
to the cutoff m^. For example, if we choose the cutoffs 
given in Fig. QJb,d), the width T(X — > J/tpcj) will be 
enlarged by a factor of 4, and the corresponding value of 
Rp/ui is smaller than 0.4. 

Our prediction of i? p / u in ([23| is a bit larger than that 
given in Ref. [19j. It is not due to a larger isospin viola- 
tion in the dispersive part of the re-scattering amplitude. 
In fact, here the isospin violation in the dispersive part 
is only about 10%, which is smaller than that expected 
in Ref. [lj| by a factor of 2. The difference is mainly 
due to the fact that the momentum factors in vertexes 
J/^D^DW and D^D^p{uj) in Fig. Q] are not consid- 
ered in the semi-quantitative estimation in Ref. [l9j |. In 
fact, the re-scattering DD* — > J/ipp(u>) is a D-wave pro- 
cess, so that the phase space smearing in (fT6|) is more 
significant than it is customarily expected. 

Furthermore, one can see from Fig. HKa) that the 
ratio Rp/DDn is roughly consistent with Eq. ((6|) ex- 
cept for the very near-threshold region where mx = 
3871.6-3871.8 MeV. However, in this region, the width of 
X — > D°D 7r°, which is obtained from Eq. (f2"TJ|) , is ques- 
tionable. Roughly speaking, the charmonium picture of 
X(3872) is not in serious contradiction with experimen- 
tal data in © if the X(3872) is slightly below the D°D*° 
threshold, i.e., mx < th n . 

In the region where mx > th n , the pure charmonium 
picture is disfavored since the prediction of R p /ddit m 
Fig. [He) is about two orders of magnitude smaller than 
the experimental data in (JfJ|). This is due to a rapid 
increase of the decay rate of X — > D a D*° + c.c. as the 
mass of X exceeds the D°D*° threshold. 

We evaluate the width by using gx = 20 GeV, and 
the result is shown in Fig. [51 One can see that T^, p = 
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35-70 KeV. Then from the obtained ratios in Fig. 4 at 
the mass slightly below th n (say around 3871.2 MeV) we 
have T(X -> rptS) = 25-100 KeV and T(X_-^ D a D°ir a ) = 
250-1000 KeV. The width of X -> D°D°j can be es- 
timated by a similar model shown in Fig. [2] and the 
value is no more than 300 KeV. Another potential de- 
cay mode of the X(3872) is the inclusive light hadron 
(LH) decay. We can use the available measurement of 
the width of IP state Xci(l-P) [6] to roughly estimate 
that T(X -> LHs) ~ T( X ci LHs) ~ 600 KeV. The 
El transition width of x cl (2P) — ► ip(2S)^ could be in 
the range 50-80 KeV; and the hadronic transition width 
of Xd(2P) -> Xci(l-P)7T7r could be 10-100 KeV. Sum- 
ming up all the estimated decay widths given above, 
we find that the total width of X(3872) is about 1-2 
MeV, which is consistent with the experimental upper 
limit Fx < 2.3 MeV. Meanwhile, the branching ratio 
Br(V — > J/ijjp) = (2-7)%, which can match the request 
of the large production rates at B-factories and at the 
Tevatron [ll, [H, HI] . 

Finally, it is worthwhile to mention that the El tran- 
sition width for \ c i(2P) — > jJ/tfj can be estimated to be 
as small as 7-11 KeV with relativistic corrections taken 
into account@, El]. Although this 2P - IS El transi- 
tion rate is sensitive to the model details duo to the node 
structure of charmonium wavcfunctions, there seems no 
difficulty in principle to explain the ratio 

IV. SUMMARY 

In summary, we re-examine the re-scattering mecha- 
nism for X(3872), as a candidate for the 2P charmonium 
state Xci(2P), decaying to J/ipp(iv) through exchanging 
D(*> mesons between intermediate states D and D* (or 
between D and D*). We evaluate the dispersive part, 
as well as the absorptive one, of the re-scattering ampli- 
tude, and find that the contribution from dispersive part 
is dominant even when X(3872) lies above the threshold 
of the neutral channel th n — m^a + TOd*o. We predict 
Rp/u ~ 1 for the mx region scanned by experiments. 
The prediction for R p /ddtt favors mx < th n and disfa- 
vors mx > thn, since in the latter case the prediction 
is two orders of magnitude smaller than the experimen- 
tal data due to a much too large decay width into real 
D°D*° mesons. Whereas when mx < th n the X can 



decay to D°D°ir only through a virtual D*° and a D°, 
and therefore the decay width of X — > D°D°ir° becomes 
much milder. Furthermore, we evaluated the width of 
X — > J/tpp with the help of a phenomenological effective 
coupling constant gx, which can be estimated from two 
different ways related to P-wave charmonium decaying 
into two charmed mesons. We find that the total width of 
the X(3872) is in the range of 1-2 MeV, and the theoreti- 
cal results for the four decay channels are roughly consis- 
tent with experimental ratios (5) and (6), as well as (4). 
The remaining problem is how to accept the low mass of 
X(3872) as the candidate of a Xci (2 ^-dominated state. 
As shown in e.g. Table I of Ref . [20j] . the mass splitting 
between Xci(2-P) and Xc2(2P) is predicted to be about 
30 MeV in potential models without the coupled channel 
effects. Including the mass shifts due to coupled channel 
effects one finds that the mass of Xci (2-P) could be further 
lowered by about 30 MeV relative to that of Xc2(2P) and 
results in a mass difference between Xc2(2P) and Xci(2P) 
of about 60 MeV (detailed discussions will be presented 
in Ref. [23]). The Z(3930) meson observed byBelle has 
been identified as the Xc2(2-P) charmonium [3J], and if 
the above estimated mass splitting makes sense, then 
the mass of Xd(2P) will be around 3872 MeV. We will 
leave the mass issue to be discussed elsewhere. Finally, 
based on the obtained results, we tend to conclude that 
a Xci(2P)-dominated state could be compatible with the 
observed decays, production in the B decay and at the 
Fermilab Tevatron, and even the mass of the X(3872). 
Therefore, aside from the molecule and other interpre- 
tations, the Xci(2-P) charmonium-dominated state could 
still be a possible assignment for the X(3872). 
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